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ABSTRACT: 

In this paper, the aeroelastic modeling and 
conceptual flutter analysis of a simple 2-D box 
wing is investigated. For each wing the torsion and 
bending elasticity is simulated by two torsional and 
longitudinal springs at its elastic axis. Furthermore, 
the connection winglet is simulated by a 
longitudinal spring. The aeroelastic governing 
equations are derived using Lagrange’s equations. 
Also, for aeroelastic loads simulation on the box 
wing, the Theodorsen aerodynamic model is used. 
The effects of various design parameters such as 
front and rear wing physical and geometrical 
properties, the connection winglet angle and 
stiffness, and aircraft altitude on the flutter of the 
system are investigated. Results show that the 
connection winglet stiffness and angle have 
significant influence on the box wing flutter 
boundary. 

1. INTRODUCTION 

The interaction between wing structural dynamics 
and aerodynamic loading is referred to as wing 
dynamic aeroelasticity. There are many dynamic 
phenomena in aeroelasticity, but the most 
important dynamic instability is flutter. Wing flutter 
is a dynamic instability of a wing associated with 
the interaction of aerodynamic, elastic, and inertial 
forces, which generally leads to a catastrophic 
structural failure of the wing [1], Because of the 
catastrophic nature of flutter, it is important that it is 
identified and avoided in the design phase. Indeed, 
the flight envelope is typically limited by the speed 
at which flutter occurs. 

Much of the immediate interest in abatement 
strategies focuses on rapid commercialization of 
green aviation technologies. Research into green 
aviation technology including reduced noise, 


emissions and fuel consumption has grown fast in 
recent years. In order to meet these goals, several 
new technology ideas and innovative aircraft 
configurations are being investigated by airplane 
designers. Innovative configurations combine the 
main parts of an aircraft (wing, fuselage, tail ...) in 
an un-conventional way with the intention to save 
on fuel burn and emissions. One such 
configuration is the Box Wing Aircraft (BWA). 

The BWA belongs to the nonplanar lifting system 
category of airplane configurations. It is a biplane 
which has a fuselage which is similar to 
conventional aircraft. Its lifting system consists of 
two wings, where the front one is swept-back and 
the rear one is swept-forward. The tips of both 
wings are connected by a vertical winglet. Figure 1 
shows the box wing aircraft model developed by 
AERO at Hamburg University of Applied Sciences. 
The most recognized benefit of this configuration is 
its low induced drag. 



Figure 1. A typical box wing aircraft [2J. 


Figure 1 shows only one possible BWA 
configuration. Many others have been proposed 
[3]. The evaluation of flutter clearance for such new 
concept wings is a major aeronautical engineering 
task. Clearly estimating the aeroelastic instabilities 
is critical to establishing the flight envelope of 
newly designed aircraft. The first recorded and 
documented case of flutter in an aircraft occurred 
in 1916 [4], After that, during the First World War, 
catastrophic failures due to aircraft flutter became 
a major design concern. The first study of the 
flutter of aircraft wings was compiled by Frazer and 
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Duncan in 1928. They introduced methods for 
analysis and prevention of flutter that are the basis 
for techniques in use today. In recent decades, 
because of the increase in wing flexibility and 
speeds, flutter analysis has become one of the 
most important aspects of wing design and, 
consequently, a lot of research has been done in 
this field. Although, the large quantity of papers 
associated with wing flutter analysis makes it 
almost impossible to consider them all, flutter 
researches on clean wings [5, 6], wings carrying 
external stores [7-9], wings with powered engines 
[10, 11], maneuvering wings [12, 13], and 

morphing wing models [14, 15] gives an idea of the 
range of research done in this area. 

Different studies on BWA designs have been 
presented in previous research. Initial 
investigations were implemented by Prandtl who 
introduced the box wing aircraft as the best wing 
system [16]. Subsequently some theoretical 
investigations have been performed by 
researchers. More recent studies have focused on 
the design of different commercial aircraft with the 
box wing configuration [17-21]. In most of this 
research some design aspects, such as aeroelastic 
constraints, were not considered in the design 
procedure. Figure 2 shows, for example, the 
recommendations plot for further improvement of 
the box wing design which was presented by 
Zohlandt [22], He presented the conceptual design 
of a high subsonic box wing aircraft and performed 
three comparison studies between box wing 
aircraft and conventional aircraft and showed that 
the box wing configuration has great potential. In 
Fig. 2, each item has been ranked on a scale from 
1 to 5 on its level of complexity and also the level 
of urgency to implement this improvement to the 
design. The resulting matrix shows the high 
necessity of aeroelastic considerations in box wing 
aircrafts design. 
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Figure 2. Recommendations plot showing urgency vs. 
complexity of the implementations needed for 
improvement of the box wing design [22], 


In recent years some works have been initiated to 
study the BWA aeroelasticity [23-26], although 
there is not yet a wide literature studying the 
aeroelasticity of Box Wings and several aspects 
need to be understood and analyzed in more 
depth. In this study, the aeroelastic modeling and 
flutter analysis of a BWA configuration is 
investigated. The results of this study can help 
designers to introduce small modifications to the 
design or come up with new BWA configurations 
that would avoid the flutter problem. 


2. AEROELASTIC MODELING 


In this research the aeroelastic conceptual analysis 
of a simple 2-D box wing as shown in Fig.3 is 
considered. In Fig. 3(a) the box wing aircraft is 
illustrated. The typical wing tip section is 
represented in Fig. 3(b), where h F and 0 F denote 
the plunge and pitch motion of the front wing and 
h A and 0 A denote the plunge and pitch motion of 

the aft wing. Wing and wingtip are considered to 
be connected with a ball joint giving no restriction 
to any relative rotational motion. 



(a) 



Figure 3. (a) Schematic of a box wing aircraft (b) 
Section A-A of the box wing. 


The aeroelastic governing equations are derived using 
Lagrange’s Equations: 
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where P and T are the strain energy and the kinetic 
energy, and Q i is the generalized force. Also, the q { are 
generalized coordinates which are defined as: 

(\-{h F 0 F h A 0 A j ( 2 ) 


The total kinetic energy of the box wing is simply: 


t=t f +t a 


(3) 


where T F and T A are the front and aft wing kinetic 
energies, respectively, which are: 
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Furthermore, the virtual work of non-conservative 
forces acting on the wing may be expressed as: 
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which leads to the following generalized forces: 
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Substituting Eqs. (3)-(8) into Eq. (1) the aeroelastic 
governing equations are obtained as: 
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The aerodynamic forces are derived from the 
Theodorsen quasi-steady aerodynamic model. Using 
this model, the aerodynamic lift and moment for the 
front wing and the aft wing can be written as: 
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Substituting Eq. (10) into Eq. (9) the following set of 
ordinary differential equations are obtained. 

Mq + Cq + Kq = 0 (11) 

Herein, M, C and K denote the mass matrix, the 
damping matrix and the stiffness matrix, respectively, 
while q is the overall vector of generalized coordinates. 


3. RESULTS AND DISCUSSION 

Following the p method, states of the box wing 
aeroelastic system are considered as: 

h F = h F exp(vt) 

6 f = 0 F exp( vl )) 

F F kv ( 12 ) 

h A = h A exp(vT) 

0 A = 0 A exp(vt)) 

Substituting Eq. (12) into Eq. (11) the following set 
of algebraic equations are obtained: 

\mA 2 +C2 + Kjq = 0 (13) 

For a nontrivial solution to exist, the determinant of 
the coefficient matrix must be set equal to zero. 
The X ’s are continuous functions of the air speed 
U. For U* 0, X is in general complex, and 
hence X = Re(/l) + / Im(A) . When Re(/L) = 0 and 
Im(2) r- 0 the wing is said to be in the critical flutter 
condition. At some point, as U increases, Re(A) 

turns from negative to positive so that the motion 
turns from asymptotically stable to unstable. 
Pertinent data for the particular wing used here are 
considered in Table 1. Also, the dimensionless 
parameter used in the numerical simulation is: 

K=k! k* (14) 
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where k* =EAH' , and EA and /' are the winglet 
tension stiffness and length, respectively. 

In Fig. 4, the modal damping of the box wing is 



plotted versus the air speed. The connection 
winglet has an angle of a = 45 deg and its stiffness 
is K = 0.5 . The flutter condition is marked by the 
crossing of the real part of one of the roots into 
positive zone. The flutter speed obtained is 
v f =81m/s, and there is modal damping in all of 

the modes below the flutter speed. 

Table 1. Characteristics of the Box Wing Model 

Parameters 


Value 

/ 


6.1 m 

V 


2 m 

b F 


0.915 m 

m F 


35.695 kg/m 

I PE 


7.172 kg m 

a F 


-0.2 

e F 


-0.1 

k hF 


1,2910 s N/m 

k 0 F 


1.6210 s N/m 

t>A 


0.732 m 

m A 


28.56 kg/m 

h A 


3.672 kg m 

a A 


-0.2 

e A 


-0.1 

k h A 


1.0310 s N/m 

k e A 


1.310 s N/m 

EA 


6.56 10 s N 

Pco 


1.224 kg/m 3 

c lf 


2 n 

c la 


In 


The flutter condition is marked by the crossing of 
the real part of one of the roots into positive zone. 
The flutter speed obtained isv^ =81m/s. It can be 

seen that there is modal damping in all of the 
modes below the flutter speed. 



Figure 4. Plot of the modal damping versus the air 
speed. 

Figure 5 shows the variation of the flutter speed 
and frequency of the box wing for selected values 


of the connection winglet stiffness due to variations 
in the front and rear wing bending stiffness. The 
connection winglet has an angle of a: = 45 deg. It 
can be seen from this figure that both flutter speed 
and frequency significantly increase by increasing 
the wing bending stiffness. Furthermore, 
decreasing the connection winglet stiffness 
improve the stability domain of the wing. 



EI(N.m 2 ) 

(a) 



EI(N.m 2 ) 

(b) 


Figure 5. Effects of the wing bending rigidity on the 
flutter boundary for selected values of the winglet 
stiffness, (a) Flutter speed, (b) Flutter frequencies. 

Figure 6 demonstrates the effect of the connection 
winglet angle on the box wing flutter boundary at 
different altitudes. The connection winglet stiffness 
is K = 0.5 . 

It can be seen that increasing the winglet angle will 
decrease the flutter speed in all altitudes. This 
means that increasing the winglet angle decreases 
the stability domain of the airplane. Furthermore, 
the effect of the aircraft altitude on the wing flutter 
speed and frequency is clearly highlighted. The 
results show an increase of the flutter speed and 
frequency with increased altitude. 

The effects of the wing length on the flutter 
speed and frequency are shown in Fig. 7. The 
connection winglet has an angle of « = 45deg 
and it is assumed that both the front and rear 
wings have the same length /. 
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Figure 6. Effects of the connection winglet angle on the 
flutter boundary for selected values of the aircraft 
altitude, (a) Flutter speed, (b) Flutter frequencies. 


(a) 



(a) 



Km) 


(b) 

Figure 7. Effects of the wing length on the flutter 
boundary for selected values of the winglet stiffness, (a) 
Flutter speed, (b) Flutter frequencies. 


As expected, by increasing the box wing length the 
flutter speed and frequency decrease. Both plunge 
and pitch effective stiffness coefficients reduce by 
increasing the wing length. 

Figure 8 reveals the effect of front to rear wing 
length ratio on the flutter speed and frequency of 
the box wing. The connection winglet stiffness is 
K = 0.5. The front and rear wing roots are 
considered fixed and the wing tips move parallel to 
the aircraft body towards the nose and tail of the 
body. The length ratio parameter is defined as: 



a (deg) 


L—If/Ia (15) 

where If is the front wing length and Ia is the aft 
wing length. Figure 8(a) shows that the length ratio 
has only a small effect on the box wing flutter 
speed. Only for high values of the connection 
winglet angle, the lower length ratios expands the 
flutter stability region of the box wing. Figure 8(b) 
shows that the front to rear wing length ratio 
influences the flutter frequency more. 

The effect of the wing chord on the flutter speed 
and frequency is shown in Fig. 9. The connection 
winglet stiffness is K = 0.5. 
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Figure 8. Effects of the connection winglet angle on the 
flutter boundary for selected values of the front wing to 
rear wing length ratios, (a) Flutter speed, (b) Flutter 
frequencies. 
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The chord ratio parameter B is defined as: 


5. RECOMMENDATIONS 


B = blb ref ( 16 ) 

where b is the wing chord and b ref is the reference 

wing chord given in Table 1. Figure 9(a) shows that 
the wing chord only has a small effect on the box 
wing flutter speed. Also, Fig. 9(b) demonstrates 
that the flutter frequency increases by decreasing 
the wing chord. 



ct (deg) 


(a) 



a (deg) 


(b) 

Figure 9. Effects of the connection winglet angle on the 
flutter boundary for selected values of the wing chord, 
(a) Flutter speed, (b) Flutter frequencies. 


4 . CONCLUSION 

Aeroelastic modeling and flutter analysis of a 
simple 2-D box wing is considered in this study. 
Plunge and pitch motions of both the rear and front 
wings are considered in the aeroelastic model. The 
connection winglet is also modeled as a 
longitudinal spring. 

The effects of various design parameters such as 
the front and rear wings length, chord, and bending 
stiffness, and also the connection winglet angle 
and stiffness, on the box wing flutter are 
investigated. The results show that increasing the 
connection winglet stiffness decreases the flutter 
speed of the box wing. Also, it is observed that 
flutter speed reduces by increasing the connection 
winglet angle. 


The 2-D approach applied here is a first step to the 
flutter analysis of the BWA. More meaningful 
results would need a 3-D approach and would 
need to use several models considering joints with 
different degrees of freedom between wing and 
winglet. According to Jemitola [27] the rigid joint 
provides the lightest structural box wing mass, 
although the rigid joint itself may be heavier than a 
ball joint. Schedl et al. [28] found the opposite: 
They obtained the lightest box wing applying a 
universal joint with unconstrained x- and z-rotation. 
The second lightest box wing was obtained 
applying a ball joint (as considered in this paper). 
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